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(57) Abstract 

A digital receiver (100) for GPS C/A-code signals. The GPS receiver (100) of the present invention provides recep- 
tion and tracking a plurality of satellites simultaneously, using four separate receiver channels. The GPS receiver (100) of 
the present invention includes an analog front-end (104) for selecting and frequency translating the received GPS signal. 
The GPS receiver (100) further includes a highspeed digital signal processor (110) for recovering the despread data of the 
GPS signal. The baseband signal is further processed by a general purpose digital signal processor (112) for signal search, 
tracking, and data recovery dperations, and a microprocessor (114) provides overall receiver control, and interface with 
the operator of the GPS receiver (100). - 
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APPARATUS FOR AND METHOD OF DOPPLER SEARCHING 
IN A DIGITAL GPS RECEIVER 

FIELD OF THE INVENTION 

Thh invention relates generally to the field of radio navigation receiver, 
OS and specifically to a ddpp,er searching «hni,ue which is adapred uTuZZ"' 
Global Posmomng S y SKm (GPS) navigation sysrem. The doppler «Lhit 

ctz; . *• **** a gps — ~ » -** «h° 1 ~ s 

10 BACKGROUND OF THE INVENTION 

Radio navigation systems are used for priding useful geographic location 
information to aircraft pilots, mariners, and even land-based vehicles such as trucks 
buses and police vehicles. Early radio navigation systems used transmitter signpost' 
^ techniques which rely on phase and timing information transmitted by several 

L5 . -terrestnal, geographically separated transmitters. One common radio navigation 
system currently in use is Loran which also relies on land-based transmitters to 
TTu ^ naVigati ° n SignaIs - latest ™«o navigation system is referred to as 
the Gtobal Positioning System (GPS), and it is maintained by the government of the 
United States of America. 

:o The GPS navigation system relies of satellites which are constantly orbicin- 

the globe. When the system is fully operational, any user of GPS, anvwhere on the 
globe„w.ll be able to derive precise navigation information including 3-dimensional 
position velocity and time of day. The system is expected to become fully 
operational in 1988 with 18 satellites in orbit. Navigation fixes using GPS are based 

5 on measurements of propagation delay times of GPS signals broadcast from the 
orbiting satellites to the user. Normally, reception of signals from 4 satellites is 
required for precise location determination in 4 dimensions (latitude, longtitude, 
altitude, and time). Once the receiver has measured the respective signal propagation 
delays, the range to each satellite is calculated by multiplying each delay by the speed 

) of light. Then the location and time are found by solving a set of four equations 
incorporating the measured ranges and the known locations of the satellites. The 
highly precise capabilities of the system are maintained by means of on-board 
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atomic clocks for each satellite and by ground tracking stations which continuously 
monitor and correct satellite clock and orbit parameters. 

Each GPS satellite transmits two direct-sequence-coded spread spectrum 
signals at L-band: a LI signal at a carrier frequency of 1.57542 GHz, and a L2 
05 signal at 1.2276 GHz. The LI signal consists of two phase-shift keyed (PSK) 
spread spectrum signals modulated in phase quadrature: the P-code signal (P for 
precise), and the OA-code signal (C/A meaning coarse/acquisition or clear/access). 
The L2 signal contains only the P-code signal. The P and C/A codes are repetitive 
pseudorandom sequences of bits (called M chips" in spread spectrum parlance) which 
10 are modulated onto the carriers. The clocklike nature of these codes is utilized by the 
receiver in making time delay measurements. The codes for each satellite are 
distinct, allowing the receiver to distinguish between signals from the various 
satellites even though they are all at the same carrier frequency. Also modulated 
onto each carrier is a 50 bit/sec data stream (also distinct for each satellite) which 
15 contains information about system status and satellite orbit parameters, which are 
needed for the navigation calculations. The P-code signals are encrypted, and 
available only to classified users. The C/A signal is available to all users* 

The operations performed in a GPS receiver are for the most part typical of 
those performed in any direct-sequence spread spectrum receiver. The spreading 
20 effect of the pseudorandom code modulation must be removed from each signal by 
multiplying by a time-aligned, locally-generated copy of it's code, in a process 
known as despreading. Since the appropriate time alignment, or code delay, is 
unlikely to be known at receiver start-up, it must be searched for during the initial 
acquisition stage. Once found, proper code time-alignment must be maintained 
25 during the "tracking" phase of receiver operation, as the user moves about. A 
mechanism for providing this alignment is called a delay-locked loop. 

Once despread, each signal simply consists of a 50 bit/sec PSK signal at 
some intermediate carrier frequency. This frequency, is somewhat uncertain due to 
the Doppler effect caused by relative movement between satellite and user, and to 
30 receiver local clock error. During initial signal acquisition this Doppler frequency 
must be searched for, since it is usually unknown prior to acquisition. Once the 
Doppler frequency is approximately determined, carrier demodulation proceeds 
using a local carrier signal derived from either a squaring or a Costas carrier 
recovery loop. In order to maintain at a constant level the dynamic characteristics of 
35 the carrier recovery and delay-locked loops as the signal strength varies, GPS 
receivers are usually provided with automatic gain control (AGC). 



05 



10 



15 



WO 87/01540 PCT/US86/01765 

-.3- . 

After carrier demodulation, data bit timing is derived by a bit synchronization loop 
and the data stream is finally detected. A navigation calculation may be undertaken 
once me signalsfrom 4 satellites have been acquired and locked onto, the necessary 
ame delay and Doppler measurements have been made, and a sufficient numberof 
data bits (enough to determine the GPS "system" time and orbit parameters) have 
been received. c 

■: To accomplish the functions described above, all known GPS receivers 

wT^r 101 teChm>IOgy a ^*~ * d ^ Passing in Z 
back-end of the receiver. The navigation computations are typically performed 

using a microprocessor, which is well-suited to the task. In addition certain 
baseband functions, such as data detection, bit timing recovery, and some Costas 
loop processing, are performed digitally in known receivers. However all 
pseudorandom code despreading, carrier demodulation, delay-locked loop 
processing, and gain control are implemented using analog components. 

An example of a prior GPS receiver is shown and described in a paper by 
Chao, Lew Cost RFILSI Technologies for Commercial GPS Receivers, Microwave 
Systems Applications Technology Conference, March 1983. Another example 
of a pnor GPS receiver is shown and described in a paper by Yiu, Crawford, and 
Eschenback, A Low-Cost GPS Receiver for Land Navigation, JOURNAL Of THE 
20 Institute Of Navigation, Fall 1982. 

There are several disadvantages to the processing approach used in prior 
GPS receivers. Due to the complex nature of the GPS signal and the complicated 
processing required, these receivers typically require large numbers of discrete 
components or highly specialized analog integrated circuits, resulting in hi-h 
25 manufacturing cost. This is especially true if the receiver is designed to process the 
four required satellite signals simultaneously, since the circuitry for one ••channel- 
must be duplicated three times. To reduce circuit complexity, some receivers 
employ what is known as sequential processing, where the hardware for one 
channel is time-shared among the four incoming signals. Receiver performance is 
30 degraded with this technique, however, since three-quarters of the information in 
each signal is lost 

Furthermore, conventional receivers suffer from the usual problems found 
with analog designs, for example, degradations in performance due to aging, 
temperature/humidity variations, and to mismatches in certain, signal path 
35 characteristics. 
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Another drawback to the processing methods employed in current GPS 
receivers is the long time needed for initial signal acquisition. As mentioned above, 
before the four satellite signals can be tracked they must be searched for in a 
two-dimensional search "space", whose dimensions are code delay and Doppler 

0 5 frequency. Typically, if there is no prior knowledge of a signal's location within 
this search space, as would be the case after a receiver "cold start", a large number 
of code delays (about 2000) and doppler frequencies (about 15) must be searched. 
Thus, for each signal, up to 30*000 locations in the search space must be examined. 
Typically these locations are examined one-at-a-time sequentially, a process which 

Iff can take up to 5 to 10 minutes. The acquisition time is further lengthened if the 
identities (i.e., codes) of the four satellites within view of the receiving antenna are 
unknown. Methods to shorten the acquisition time have been devised but are quite 
expensive to implement. One technique, for example, employs surface acoustic 
wave filters matched to each of the 18 satellite codes to effectively perform the 

15. despreading. Another technique utilizes multiple conventional despreading circuits 
operating in parallel in order to search several code delays simultaneously. 

It is apparent from the above discussion that conventional GPS receivers, 
which rely mainly on analog technology, suffer from many disadvantages. Most 
prior receivers perform the GPS signal search in a sequential manner. The usual 

20 procedure is for the receiver to make a guess as to the correct code delay, setting the 
C/A code generator accordingly. Then, for that code delay, all frequency bins are 
examined for presence of signal. If no signal is found. , a code generator is 
incremented to the next code delay and the process is repeated. The detection 
method typically employed in checking for signal presence can be loosely described 

25 ; as a variables-frequency average power "meter" followed by a threshold decision. 
For a given code delay /Doppler frequency trial, the C/A signal is multiplied by the 
locally generated C/A code, and mixed with a sinusoid whose frequency is equal to 
the current Doppler guess. The resultant signal is filtered by a filter whose 
bandwidth is commensurate with the frequency search interval Af. If the desired 

30 signal is present in the bin being searched, the output of said filter will contain a 
despread 50 bit/sec data signal. The output of the filter is squared and averaged to 
obtain the average power, which is then compared to a threshold. If the average 
power is greater than the threshold, the signal is deemed present in the current bin. 
Otherwise, the signal is deemed absent, and the search proceeds to the next bin. A 
3 5 search technique which typifies a conventional practice is described in a paper by 
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Spilker, J.J., Jr., Global Positioning System: GPS Signal Structure and 
Performance Characteristics,. THE INSTITUTE OF NAVIGATION. Vol. I, 'pp. 29-54 
1 980* ~ * * ? 
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. Due to the large number of bins to be searched (around 30 000) 
conventional receivers require many minutes to find a signal if a sequential search 
method is employed. Additional time may be incurred if the identities (i e codeo 
of the satellites within view of the receiving antenna are unknown. Methods to 
shorten the search time have been devised but are quite expensive to implement. 
Generally, these are based on the idea of examining several code delays 
szmultaneously for each Doppler frequency in order to speed up the process For 
example, one technique employs multiple despreading circuits, each with its own 
code generator and average power detector, operating in parallel. Therefore, a need 
exists for a doppler search technique which allows the GTS receiver to rapidly and 
efficiently search availble doppler frequencys and to improve the lock-up tiame in the 
15: GPS receiver. 

SUMMARY AND OBJECTS OF THE INVENTION 

In summary. The present invention contemplates a GPS search technique 
for the reception of GPS C/A-code signals from up to four satellites simultaneously 

20 using four separate receive channels. The principles of the present invention can also 
be applied to simultaneous reception of GPS P-code signals; to sequential reception 
of C/A or P-code signals using one time-shared channel: and to reception of 
direct-sequence-coded spread spectrum signals in general. 

In accordance with the preferred embodiment of the invention, the search 

25 technique of the present invention performs doppler search in a spread-spectrum 
receiver, for receiving one or more direct-sequence-coded signals, each of which 
may be Doppler shifted from its nominal carrier frequency, wherein for each signal, 
a code delay/Doppler frequency search procedure must be carried out prior to 
activation of delay-tracking and carrier-recovery loops. For each code delay search 

30 trial, the present invention performs the steps of inputting N consecutive digitized 
samples of a signal equal to the filtered product of the down-converted received 
signal multiplied by the locally-generated code of the desired signal, the code delay 
of the locally-generated code being that of the current code delay search trial, the 
filtering passing, the band of frequencies corresponding to the entire range of 

3 5 possible Doppler-shifted carrier frequencies of the down-convened received signal: 
The method of the present invention then calculates the N-point discrete Fourier 
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transfonn of the block of N samples to produce N complex-valued transformed 
output samples, each corresponding to one of N equally-spaced search frequencies 
in the band from zero Hertz to fs (or, equivalentiy, in the band from -f s /2 to +f s /2), 
f s being the sampling frequency of the samples input in the previous step. The 

05 square of the magnitude of each of said N complex-valued output samples is then 
calculted to produce N power quantities, corresponding to the N said search 
frequencies. Steps a through c are performed M times, summing the M power 
quantities produced for each of the N said search frequencies, to produce N 
quantities proportional to average power. The maximum of the N average power 

la quantities is noted and the corresponding frequency is calculated. The maximum 
average power quantity is compared to a threshold, the fact of said threshold being 
exceeded signifies that the current code delay is approximately correct (local code 
approximately time-aligned to received code), with the frequency of maximum 
average power found in step e indicating the approximate Doppler frequency shift, 

IS and the fact of said threshold not being exceeded signifying that the current code 
delay is not approximately correct, and that another code delay should be tried. 

Accordingly, it is one object of the present invention to provide a GPS 
search technique which rapidly locates a GPS signal. 

It is another object of the invention to provide a GPS search technique 

20 where all possible doppler signals are searched simultaneously. 

It is yet another object of the present invention to provide an improved 
means for rapid initial acquisition of GPS signals. 

It is still another object of the present invention to provide a GPS search 
technique which can be readily imptmented in a general purpose Digital Signal 

25 Processor. . 

These and other object^ of the invention will become apparent to those 
skilled in the art upon consideration of the following description of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
3Q Figure I is a block diagram showing the architecture of the GPS receiver of 

the present invention. 

Figure 2 is a detailed block diagram of the analog front end 104 of Figure 1, 
and showing connections to the antenna and A/D converter. 

Figure 3 is a detailed block diagram of the first DSP section 1 10 of the GPS 
35 receiver of Figure 1. 
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Figtues 4A and 4B are schematic diagrams of alternate embodiments of quadrature 
nuxer/decimator310ofFigure3. quadrature 

Figures 5A and 5B are schematic diagrams of alternate embodiments of 
Iowpassfdter/decimator322 ofFigure3. wiments of 

106. FiSUW ^ " ' SPCCtraI Pl0t ° f * ° A SignaI " *° ° UtpUt ° f c °nvener 

Figure 6B is a spectral plot ofaC/A signal at the output of quadrature 
mixer/decunator 310. Hu*orature 

Hgure 6C « a spectral plot of a despreadda^ signal, and 
the responses of the filters of Figures 5 A and 5B. hereto. 

Figure i Fi8Ure 7 " a det ^ led bIock of the second DSP section 112 of 

Figure 8 is a detailed block diagram of the operations performed by 
tracking/data processor 704 of Figure 7. 

Figures 9A, 9B, and 9C illustrate the problem of searching for a GPS C/A 
signal over many code delays and Doppler frequencies. 

Figure 10 is flow diagram of the FFT-aided maximum spectral power 
measurement performed by search processor 702 of Figure 7 in accordance with the 
method of the invention. 

Figure 1 1 is a flow diagram of the search procedure performed by 
microprocessor 114 of Figure 1 in conjunction with search processor 702 in 
accordance with the improved search method of the invention. 

DETAILED DESCRIPTION OF THE DRAWINGS 

Figure I is a block diagram showing the architecture of the digital GPS 
receiver of the present invention. The receiver 100 includes an antenna 102 coupled 
to an analog front-end 104, for the purpose of receiving and selecting the desired 
composite GPS signal (comprising signals from several satellites) at the LI 
frequency and for modifying the signal such that it is suitable for conversion to 
digital (binary) form. This signal modification includes filtering, ampUfication. and 
frequency translation. The analog front-end 104 is coupled to A/D converter 106 
which samples and converts the received analog signal to the digital form required 
for subsequent receiver processing. An oscillator 108, coupled to analog front-end 
104 and to A/D converter 106, supplies a clock signal at an appropriate sampling 
frequency. In the context of the present invention, this clock frequency is to be 
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carefully chosen, as will be discussed further below. In the preferred embodiment 
of the invention, the clock frequency is approximately 38. 192 MHz. 

The output of A/D converter 106 is coupled to a first digital signal processor 
(DSP) 110. The DSP 110 provides quadrature mixing to approximately zero 
0 5 frequency, C/A code generation and despreading, early-late processing for deriving 
auxiliary signals for delay-locking, lowpass filtering, and sample rate reduction. 
According to the principles of the present invention, the DSP 1 10 processes received 
signals from 4 satellites simultaneously. Oscillator 108 is also coupled to DSP 110, 
and supplies the clock signal at the fundamental sampling frequency. First DSP 1 10 
la operates at a relatively high speed, owing to the high sample rate of the digitized 
signal supplied by A/D converter 106. However, in spite of the high operating 
speed, and in spite of the large amount of processing required, the operations of 
DSP 1 10 may be simplified, in accordance with the principles of the present 
invention, to such an extent that DSP 1 10 is readily adaptable to implementation in a 
L5 single monolithic integrated circuit (IQ. 

First DSP 1 10 is coupled to a second DSP 1 12 which performs signal 
search, tracking, and data recovery operations. The DSP 1 12 also performs signal 
processing for four channels simultaneously. These operations include Costas loop 
carrier recovery and demodulation, delay iocked-ioop filtering and CI A code delay 
20 control, data bit synchronization and data detection, and Mr l -aided spectral power 
* measurement for search. Owing to the relatively slow sample rate of the signals 
supplied by first DSP 110 (by virtue of the lowpass filtering and sample rate 
reduction performed therein), the operations of second DSP 112 are readily 
implemented using general-purpose, programmable DSP ICS which are readily 
25 ' available from a number of manufacturers. 

Second DSP 1 12 is coupled to a microprocessor ftiP) or \x? system 1 14. 
which handles overall receiver control and navigation functions. These functions 
include interpretation of data messages and other parameters measured by DSP 112 
(e.g., code delays and Doppler frequencies), navigation fix computations, and 
30 selection of the four satellites to be tracked. Microprocessor 1 14 also performs, in 
conjunction with second DSP 112, an improved search algorithm to be described in 
more detail below. The functions of pP 1 14, with the exception of the improved 
search algorithm of the present invention, are typical of those performed by other 
GPS receiver back-end \xP\ and will not be further elaborated. An example of the 
35 operation of the back-end microprocessor 1 14 is disclosed in an article by Borei M. 
L, et al., Texas Instruments Phase I GPS User Equipment, THE INSTITUTE OF 
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NAVIGATION: THE GLOBAL POSITIONING SYSTEM, pp. 87-102, 1980. Details of a 
navigation fix algorithm in particular may be found in an article by Noe P S and 
Myer, K. A., A Position Fixing Algorithm for the Low-cost GPS Receiver Ip Pp 

O, ^<^° NS ° N AEROSPACE **> ^CTRONIC SYSTEMS, AES- 12 Vol 2 
OS pp 295-297, March 1976. ' oL 2 " 

Figure 2 shows a more detailed block diagram of analog front-end 104 
Analog front-end 104 includes a preselector 202 coupled to antenna 102, for the 
purpose of selecting and amplifying the desired composite GPS signal Signal 
selection is accomplished by the preselector which incorporates a bandpass filter 
10 appxox^ately centered at the LI frequency of 1575.42 MHz. The parameters of the 
filter (bandwidth and number of poles) are chosen as a compromise between the 
confl ctmg requirements of good interference rejection and low signal distortion, as 

If ops'™ ^ " " t In * Presem invcntion ' which is directed » 

L5 Those skilled m the art will appreciate that ceramic filters are particularly well suited 
for use as a preselector in the context of the present invention. 

^ ^ ° UtPUt ° f PrCSeleCt0r 202 is cou P Ied "> one input of an analog mixer 
204. The second input of analog mixer 204 is coupled to the output of harmonic 
generator 210, which supplies a substantially sinusoidal signal exhibiting a 

20 frequency equal to N times the sampling frequency f s . The input of harmonic 
generator 210 is coupled to oscillator 108, which supplies a clock signal at 
frequency f s , whose harmonic is to be generated. Harmonic generator 2 10 may be 
reahzed using any of several standard techniques, for example, a step recovery diode 
may be used as a pulse generator, which is followed by bandpass filtering to pick 

25 off the desired (Nth) harmonic. " 

According to the method of the invention, the sampling rate f s and the 
harmonic number N are chosen so as to satisfy the following relationship: 

f Ll = 1575.42 MHz = (N .25) • f_ 



30 



If this is satisfied, then the difference-frequency signal appearing at the output of 
mixer 204 will be approximately centered at the frequency f s /4, which will be 
explained in more detail below. Thus, the combination of harmonic generator 210 
and mixer 204 act to frequency-translate' the incoming signal from a carrier 
35 frequency of 1575.42 MHz to a carrier frequency equal to approximately f s /4. In the 
preferred embodiment of the present invention, N=41, f s =38.192 MHz, and the 
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plus (+) sign in the above equation is used In this case, the GPS signal is translated 
to approximately 9.548 Mhz. 

The output of mixer 204 is coupled to the analog input port of an 
analog-to-digital (A/D) converter 106. Oscillator 108, also coupled to an A/D 
05 converter 106, supplies a clock signal which causes the A/D to sample, and digitize 
the input waveform at the sampling rate f s . The multWrit digital output of A/D 
converter 106 is coupled to DSP 1 10 for further processing. In order to avoid the 
need for a separate sample-and-hold circuit, A/D converter 106 is preferably of the 
"flash" variety. Several flash A/D converters are available off-the-shelf from various 
10 manufacturers. For example, the TDG1029 6-bit A/D converter available from TRW 
may be utilized. TRW LSI Products Division, TRW Elecronic Components Group, 
P.O. Box 2472, La Jolla, CA 92038. 

The quantization noise introduced by A/D converter 106 adds to the antenna 
and other front-end noises produced by preselector 202 and mixer 204, It is 
IS desirable to adjust the gain of preselector 202 so that this quantization noise does not 
significantly degrade the overall front-end noise figure. The resolution (number of 
bits) of A/D converter 106 should be chosen to provide adequate dynamic range for 
the types of signals that will be received, including desired GPS signals and 
interferers. In the preferred embodiment of the present invention, a 6-bit converter 
20 is utilized. The sampling rate f s must be large enough to provide adequate protection 
against aliasing, in accordance with the Nyquist sampling theorem. As stated 
previously, a sampling rate of 38.192 MHz satisfies this criterion, as well as the 
relational constraint described above. 

Note that analog front-end 104 represents a very minimal amount of circuitry. 
25 in contrast to the analog sections of prior receivers. This has significant beneficial 
implications with regard to manufacturability, and immunity of receiver performance 
to environmental variations. The output of A/D converter 106 is coupled to first 
DSP 110, which is depicted in Figure 3. DSP 110 includes a quadrature 
mixer/decimator 3 10, coupled to A/D converter 106, which accepts as input the 
30 sampled, digitized composite GPS signal comprising C/A signals from nominally 
four satellites, plus front-end noise. Mixer/decimator 310 performs quadrature 
mixing to translate the composite signal by f s /4, from a center frequency of 
approximately fs/4, down to a carrier frequency of approximately zero Hz. 
Although translation to zero Hz is not strictly necessary, it is preferred since it 
35 simplifies subsequent processing. After this frequency translation, mixer/decimator 
310 performs very simple lowpass filtering, followed by sampling rate reduction, or 
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"decimation" The quadrature mixer/decimator 310 will be discussed in more detail 
below in conjunction with Figures 1 and 2. 

As is well known in the art, quadrature mixing produces com»u* 
number ^ outputs with two components: an in-phase (I) Lponen. and L 
out-of-phase or quadrature (Q) component The I and Q outputs Q f quadrati, 
mixer/decimator310are eachcouol e df rt f rt „, u , quadrature 
360 a „H «n ! coupled to four channel processing units 320, 340 

360, and 380 corresponding to channels 1, 2, 3, and 4, respectively The/' 
processing units perform code despreading and related tasks for four'salmte 
stgnals. Processmg units 320, 340. 360,. and 380 are identical in stn.ctuL- hen 1 
only the workings of channel 1 processor 320 wiU be described. 

Channel processor 320 includes a C/A code generator 330, which produces 
two signa l) a "prompt" signal, which is a duplicate of the C/A codetf 
deseed satemte, and 2) a "late-early" signal, which is defined as one-half the 
anthmetxc difference of the prompt signal delayed (late) by a time interval T d , and 
the prompt signal advanced (early) by T d . In mathematical terms, if the prompt 
signal is denoted by P(t), then the late-early signal is: 

LE(t) = { P(t-T d ) - P(t+T d ) } / 2 

. Since the prompt signal is a C/A code and hence takes on the values .1 and 
-1, the late-early signal can be seen to take on the three values .1, -1, and 0 Code 
generator 330 can be designed according to standard techniques- C/A code 
generation is fully described, for example, in Air Force document SS-GPS-300B 
System Specification for the VAVSTAR Global Positioning Svstem. The required 
signals may be obtained by passing the basic C/A code through a shift renter 
which is tapped at the input, center, and output to provide the earlv, prompt, and late 
signals, respectively. The time interval T d is chosen to yield the desired 
delay-discriminator characteristic, or "S-curve", as is well-known in the art 
Typically, T d is approximately one-half the duration of a chip, or in this case about 
0.5 microseconds, since the C/A code chip rate is 1.023 Mchips/sec. The control 
input of C/A code generator 330 is coupled to second DSP 1 12, which transmits 
instructions as to which of the 18 possible C/A codes is to be generated, and as to 
the desired timing relationship (code delay) between the generated code and the code 
of the received satellite signal. 

The prompt output of C/A code generator 330 is coupled to one input each 
of multipliers 332 and 334. The second inputs of multipliers 332 and 334 are 
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coupled to the I and Q outputs, respectively, of quadrature mixer/decimator 310. 
Assuming proper time-alignment, multipliers 332 and 334 effectively perform 
despreading, or removal of the C/A code, from the signal whose code matches that 
produced by C/A code generator 330. Thus, assuming proper time-alignment, the 
05 I-Q signal appearing at the output of multipliers 332, 334 is a50 bit/sec data signal. 
The carrier frequency of this data signal is typically not exacdy zero Hz, rather, it 
varies from zero Hz according to the Doppler shift caused by satellite and user 
movement, and by drift of oscillator 108 from its nominal frequency. For typical 
land users, the Doppler shift due to satellite and user movement amounts to about +/- 
10. 2.7 parts-per-million (ppm). Assuming an oscillator drift of +/- 2.0 ppm, the total 
Doppler shift is typically +/- 4.7 ppm, or about +/- 7.5 kHz with the nominal LI 
carrier frequency of 1575.42 MHz. Since the I-Q signal at' the output of multipliers 
332, 334 is known to lie within the frequency range -7.5 kHz to +7.5 kHz, it may 
be lowpass filtered to remove unwanted noise outside that band. In turn, the 

15: sampling rate of the signal may be reduced in accordance with Nyquist's sampling 
theorem. These are the functions of lowpass filter/decimators 322 and 324, whose 
inputs-are coupled to the outputs of multipliers 332 and 334, respectively. The PI 
and PQ (prompt I and Q) signals, output from lowpass filter/decimators 322 and 
324, respectively, are coupled to DSP 112 for further processing. 

ZO. Figures 6a, 6b, and 6c show the magnitude spectra of a GPS OA signal as 

it is processed through DSP 1 10. Figure 6a shows the one-sided spectrum of a C/A 
signal input from A/D converter 106. It exhibits a (sin x/x) type response with the 
first nulls located approximately 1 MHz from the center frequency of fs/4 (9.548 
Mhz in the preferred embodiment). This shape is due to the C/A code, whose chip 

25 rate is 1.023 MHz. Figure 6b shows the spectrum of the translated I-Q signal 
appearing at the output of quadrature mixer/decimator 3 10. The spectrum of an I-Q 
signal is formed, in general, by taking the Fourier transform of a complex signal 
whose real part is the I signal, and whose imaginary part is the Q signal. Note that 
the translated C/A signal is centered at approximately zero carrier frequency. 

30 Figure 6c shows the spectrum of the prompt I-Q signal after despreading. 

Note that the spectral spreading caused by C/A code modulation has been removed, 
leaving the 50 bit/sec data signal centered at the Doppler offset frequency. 

Referring again to Figure 3, the late-early output of C/A code generator 330 
is coupled to one input each of multipliers 336 and 338. The second inputs of 

35 multipliers 336 and 338 are coupled to the I and Q outputs, respectively, of 
quadrature mixer/decimator 3 1 0. The I-Q signal appearing at the output of 
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multipliers 336, 338 is an auxiliary -differential" signal which is utilized later in the 

the afferent*! I-Q signal is confined in frequency to the range of\Lible £ , 
offsets, that is, to the band between -7.5 kHz and + 7iiST£*[ 
0 5 owpass filter/decin.tors 326 and 328, coupled to the outputs n^££ 
338, respectively, remove unwanted energy outside this band, and reduce *e 
samp mg rate accordingly. The pi and DQ (differential I and Q) signal ou»u! 
ron^wpassmter^ec^ 



As mentioned previously, channel processors 340, 360, and 380 are 
identical « structure to channel processor 320. During normal operation each 
processing unit separates, from the received 

composite signal, the particular satellite signal whose C/A code matches that 
15 produced by that unit's code generator. 

A clock/control generator 390 is provided as part of DSP 1 10 for generating 
clock and control signals required by the various digital processing, circuitries 
Oscillator 108, coupled to clock/control generator 390, provides the basic clock 
sxgnalfrom which the other clock and control signals are derived. The architecture 

IncJ f -'^ ° f * C PreSem inVCnti0n HaS bee " *° that the circuitry 

of DSP 110 is amenable to realization on a single mo,.. chic integrated circuit This 
is chiefly due to the fact that complicated, multi-bit di,uol multiply operations, which 
are typically necessary in digital signal processing, are avoided throughout. As is 
well-known in the an, multi-bit multiplications (e.g.. multiplication of two 8-bit 
25 numbers) are quite costly in terms of circuit area, power supply current 
consumption, and execution time. This is especially true at high operating speeds or 
samphng rates, such as are required in a GPS receiver. Multi-bit multiplications are 
avoided in DSP 110 by several methods. First, recall that the prompt signal supplied 
to multipliers 332 and 334 exhibits only the values +1 and -1, and that the late-earlv 
30 signal supplied to multipliers 336 and 338 exhibits only the values +1, -1, and 0. 
Those skilled in the art will appreciate the ease with which multiplication by such 
simple factors may be implemented. Only a few logic gates are needed per bit of 
word length of the I and Q signals being multiplied. The OA code generator 330 
contains no multipliers, and may be easily configured using gates and flip-flops, as 
35 shown in the document cited above. Multi-bit multipliers may be eliminated from 
lowpass filter/decimators 322, 324, 326, and 328 by the use of special filter 
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structures, as will be explained later Finally, such multipliers are avoided* in 
quadrature mixer/decimator 310 by virtue of several architectural features of the 
receiver of the present invention, including a special choice of quadrature mixing 
frequency, which will now be discussed. 

0 Si Recall that the main function of quadrature mixer/decimator 3 10 is to 

translate the composite signal input from the A/D converter down to zero frequency. 
In.general, quadrature mixing is performed by multiplying the signal to be translated 
by cos Wgt and siii Wgt, where Wg is the desired amount of radian frequency 
translation. Mathematically, this may be thought of as multiplying by the complex 

Id quantity exp{jWgt} = cos Wgt + j sin W Q t, and regarding the real and imaginary 
parts of the product as the I and Q outputs, respectively. In a digital implementation, 
the cosine and sine waveforms are sampled, i.e., the time variable t becomes the 
discrete time variable nT r where T = l/f s * 1/sampling rate, and n is a sample count 
integer (1,2,3,4, ...etc.) With regard to circuitry, the digital implementation requires 

15 some means to generate the cosine and sine values, such as a ROM table, plus two 
multiplier circuits. In general, for good performance the cosine and sine values must 
be accurately represented in a substantial number of bits, which means that the 
cosine/sine multipliers are necessarily of the complicated multi-bit variety. 

Multi-bit multipliers and cosine/sine generation circuitry are avoided in the 

20 quadrature mixer of the present invention by a special choice of translation 
frequency. Recall that the composite signal input from A/D converter 106 is 
centered approximately at the frequency f s /4, by virtue of the previously-described 
constraint specifying the relation between the sampling frequency f s and the 
harmonic generator parameter N. Thus, to translate the composite signal to zero 

25c frequency from f $ /4, the translation frequency Wg must be 2;tf s /4. This means that 
the sampled complex mixing waveform exp{W 0 t} becomes exp{jn7T/2}. which is 
simply of the form ...1, j, -1, -j, 1, j, -1, -j, ...etc. The sampled cosine and sine 
waveforms are thus simplified to ...1, 0, -1, 0, and ...0, 1, 0, -I, respectively. 
Those skilled in the art will appreciate the ease of generation of such simple 

30 sequences, and of multiplication therewith. In particular, the necessity for 
complicated multi-bit multiplier circuits is eliminated. 

The above-described simplification of the quadrature mixing operation is 
made possible fundamentally by an architectural feature of the GPS receiver of the 
present invention. All known GPS receivers perform carrier demodulation, or 

35 Doppler removal, either prior to, or in association with, despreading. Such an 
approach, if adopted in a digital receiver, would imply quadrature mixing to exacdy 



WO 87/01540 . PCT/USS6/01765 

-15- 

aero Hz, rather than to approximately zero Hz, as is the case with the „„, 
invention. That is, the translation frequency would be, instead of f ,4 f T 
f Doppler Sincethe * • r s'4 + 

Doppler frequency is generally different for eaeh satellite signal, , his would 

05 nee. S s.u« four separate quadrat ntixing operations. More imp~ *e 
spectal retadonshtp between the nans.ation frequency and sampUng , at e woual 
destroyed, e.unsnating the possibility of simpUfyiag Ae ^ LefZL ana 
avo.d.ng ra ^. b i tm „ lti p Uercin:uits . By con Mst , the receiver of ttTesent 

10 Idl ^"^^'^"^^Pteading. in fact, aftertheS 

10 ^^^-oP=ra nOT swh i chfoUowdespr^ i „ g .Th iS panitio„ingoffu„cr„s 
a^owscatxterdentodulation to be perforated in DSP 1 12 using ptograLtab.e DSP 
K whtch are wel.-suited to the tasks of waveform generation and multi-bi, 
multiplication. It also makes possible the above-mentioned simplifications, allowing 
quadrature moxer/decimator 310 to be easily implemented, as will now be described. 

Ftgures 4a and 4b show two alternate embodyments of quadrature 
mtxer/decmator 3 10. Both embodiments perform quadrature mixing by the 
comp.ex sequence ...1, j, -1, -j,,.. etc., followed by simp.e filtering and decimation 
operanons. The filtering is performed ,o remove the negative-frequencv image of 
the composite GPS signal, which is centered a. the frequency ui after the 
20. quadrature mixing process. Sincethe signal bandwidth is reduced by such filtering 
samplmg rate reduction (decimation) is also performed to allow subsequent circuitry' 
to operate at lower speed. The circuit of Figure 4a effectively performs quadrature 
mixing followed by a 2-tap digital lowpass filter (z-transform transfer func!ion H(z) 
- 1 *;z ) and decimation by 2. The circuii of Figure 4b in effect performs 
2= quadrature ntixing followed by a 4-tap lowpass filter (transfer function Hfz) - 1 - 
z + z . +• z* J ) and decimation by 4. 

The embodiment of quadrature mixer/decimator 3 10 of Figure 4a includes two 
multi-bit flip-flops or latches 402 and 412 (denoted F/F), each coupled to A/D 
converter 106, for the purpose of latching the multi-bit signal samples emanating 
30 from the A/D and separating the sample sequence into two alternating half-rate 
sequences. Flip-flops 402 and 412 are both provided clock signals of frequencv 
f s /2. However, the clock signal for flip-flop.412 is 180 degrees out-of-phase with 
that for fli P -fl op 402, so that while flip-flop 402 latches even-numbered input 
samples, flip-flop 412 latches odd-numbered samples (or vice versa). The clock 
35 phase inversion may be 
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conveniently implemented: by, for example, inverter 411. Flip-flop 4 14, coupled to 
the output of flip-flop 412 and provided the same clock signal as flip-flop 402, acts 
to time-align the two sample streams. Multipliers 406 and 416, coupled to the 
outputs of fliprflops 402 and 414, respectively, mix each stream with alternating 

05 sequences. Multipliers 406 and 416 can be easily implemented due to the 

simple nature of the factors +1 and -1. The outputs of multipliers 406 and 416 
represent the I and Q outputs, respectively, of quadrature mixer/decimator 310, and 
are coupled to channel processing units 320, 340, 360, and 380 for further 
processing. The sample rates of the I and Q outputs are reduced by half from that at 

10 the input, down to f s /2 (19.096 MHz in one receiver embodiment). 

The embodiment of quadrature mixer/decimator 3 10 of Figure 4b includes 
flip-flops 442, 452, and 454, and inverter 451, connected in a manner identical to 
that of circuits 402, 412, 414, and 411 of Figure 4a, and producing the same effect. 
That is/circuits 442, 452, 454, and 451 serve to latch the sampled signal coming 

iq- from A/D 106, and to separate it into two synchronized, half-rate sample streams. 
One sample stream, appearing at the output of flip-flop 442, is coupled to flip-flop 
446. Flip-flop 446 is provided the same clock as is flip-flop 442 (at frequency f $ /2), 
and provides a delay of one sample period. Subtractor 448 subtracts the delayed 
samples appearing at the output of flip-flop 446 from the undelayed samples 

20 appearing at the output of flip-flop 442. The output of subtractor 448 is latched by 
flip-flop 450 at the rate f s /4, as instructed by a clock signal at frequency f s /4. Every 
second output of subtractor 448 is thus discarded. Identical processing is performed 
on the sample stream output from flip-flop 454, by flip-flops 456 and 460 and 
subtractor 458. The outputs of flip-flops 450 and 460 represent the I and Q outputs, 

25 respectively, of quadrature mixer/decimator .310, and are coupled to channel 
processing units 320, 340, 360, and 380 for further processing. The sample rates 
of the I and Q outputs are reduced by one quarter from that at the input, down to f s /4 
i (9.548 MHz in one receiver embodiment). Note that this rate is one-half that 
provided by the circuit of Figure 4a, allowing subsequent circuitry to operate at a 

30 lower rate. Although somewhat more complex, the circuit of Figure 4b would be 
preferred over that of Figure 4a for realizations where processing speed or power 
consumption is critical. 

The foregoing discussion describes how the operations of quadrature 
mixer/decimator 310 may be substantially simplified, allowing for easy 

35 implementation. Attention is now turned to the post-despreading iowpass 
filter/decimators 322, 324, 326, and 328. Digital filters in general require three 
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types of components: 1, Iatt hes or ffiplnpps, which provide delay, 2) adder, or to 

WbaXK>n ' Md 3) MUtfpita, are used to ~ 2 

signal samples flowing through a digital film, by various gain consl " 
eon™, the fdte* nansfer function, or frequency response. The 
comptecty depends on the signal sample and gain coefficient word .engo* Z 

Zr ^ ' h : d r niC ^ ** i" ^rnuncTof * 

filler, The coefficnt word leng,h affects the ihe degree to which the fitted 

frequency response can be made to approximate some desired response. C 
coefflcten. quannaadon (i.e.. short coefficient word length, is des^e £Z 
tmpfcrnentano, standpoint but can cause substantia! distortion of the respoTe n 
some filter structures. This is especiaUy nue if me passband of the respond ve£ 
narrow compared to me satnphng frequency, as is the case here. To overcome t„7s 
pmMem, the mcetverof the present invention emptoys special fflter structures which 
allow me coefflcents to be very coarsety quantized, to such an extent that muhiplier 
15- circuits are ebrninated entirely. 

Figures 5a and 5b detail two alternate embodiments of multiplierless digital 
filters which may be used as the lowpass filter/decimator 322 (and also 3?4 326 
and 328, smce they are identical). Circuit 5 10 of Figure 5a can. basically be' 
described as an M-sample sum-and-dump filter. As is well known in the art' the 
function of such a filter (also called an integrate-and-dump) is to, for every set of M 

consecutive input samples, form the sum of the samples and output the result This 
is equivalent to a finite impulse response (FIR) filter with M equal gain coefficients 
whose output is sampled once for every M input samples (decimation by M) This' 
function may be conveniently implemented by circuit 510, which comprises an 
accumulator formed by binary adder 502 and latch or multi-bit flip-flop 504. and an 
output latch 506. Input samples, arriving at sampling rate f r . are summed bv adder 
502 into accumulator latch 504, which is clocked at the same rate f r . After M input 
samples have been summed, the result of that sum, appearing at the output of latch 
504, is saved by (or "dumped" into) output latch 506, which is provided a clock at 
the decimated sampling rate tyM. The output of latch 506 is the lowpass filtered 
output signal whose sampling rate has been desirably reduced by the factor M to 
fr/M. Accumulator latch 504 is also provided with a clear, or reset, signal at the rate 
fl/M, which causes it to be cleared to an all-zero state after each M-sample sum has 
been formed and output to latch 506. the accumulator is thus readied for 
summation of the next M input samples. In the context of the present invention, the 
input sampling rate f r is equal to either f s /2 or f $ /4, depending on whether the circuit 
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of Figure 4a, or that of Figure 4b, is utilized for quadrature mixer/decimator 310. 
Thus, if f s is 38.192 MHz, fj is either 19.096 MHz or 9.548 MHz. 

The magnitude vs. frequency response of sum-and-dump 510 exhibits a sin 
(x/x) shape centered at zero frequency, an example of which is shown in Figure 6c. 
05 For a given input sampling rate fj, the bandwidth of the filter is controlled by the 
parameter M, with the first stopband null appearing at the frequency fj/M. The 
bandwidth should be chosen as a compromise between the conflicting requirements 
fori) low signal-to-noise ratio (SNR) degradation for despread data signals within 
the Doppler range +/- 7.5 kHz (implying large bandwidth), and 2) low output 
L0: sampling rate, so that the processing burden on DSP 1 12 is minimized (implying 
narrow bandwidth). 

Although simple to implement (no multipliers), a sum-and-dump filter 
possesses somewhat undesirable response characteristics, namely, a rather rounded 
passband, and a stopband which rolls off rather slowly with frequency. These 
15 characteristics make for somewhat degraded performance compared to that of more 
elaborate filters. The ideal filter for this application would be a rectangular or 
"brick-walT filter with a cutoff frequency equal to the maximum expected Doppler 
shift of about 7.5 kHz. Although not practicality realizable, such a filter would 
allow the output sampling rate to be reduced down to about 15 kHz, without causing 
20 the SNR degradation normally produced by noise aliasing and signal attenuation. 
By contrast, a sum-and-dump filter designed for the same output sampling rate 
introduces SNR degradations of up to 3.9 dB, with the worst degradation occuring 
for signals at the Doppler limits. For most GPS receiver applications this level of 
degradation is intolerable. Therefore, in most cases sum-and-dump 510 would be 
25 designed with a wider bandwidth, the penalty, of course, being a higher output 
sampling rate and a larger processing burden imposed on DSP 1 12. 

In cases where the SNR/output sampling rate tradeoff allowed by 
sum-and-dump 5 10 is insufficient, an alternate filter embodiment of the type shown 
in Figure 5b may be employed. The filter of Figure 5b exhibits a frequency 
30 response more closely approaching that of an ideal filter, yet is easily implemented 
as it contains no multiplier circuits. Briefly, iowpass filter/decimator 322 of Figure 
Sbcan be described as the cascade of an Ml-sample sum-and-dump, a second order 
multiplieriess recursive filter with two poles and two zeroes, and an M2-sample 
sum-and-dump. The overall filter may be viewed as an MlM2-sample 
3 5 sum-and-dump which has been split into two sections, between which a 2nd order 
recursive section is placed which acts to "square up" the passband response and 
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improve me stopband attenuation. Sampling rate redaction, or decimation is 
accomplished in two stages, first, by the factor Ml by the first sum-and-dump and 
agam by the factor M2 by the second sum-and-dump, yielding an overall decimation 
factor of Ml x M2. Thistwo stage decimation allows the second and third filter 
05 sections to operate at a reduced clock rate, thereby reducing power consumption 

In detail, Iowpass filter/decimator 322 of Figure 5b includes a first 

(either f /2 or f s /4, as discussed above), and performs sum-and-dump operations 
. identical to those described above, except that the parameter M is now Ml. Filtered 
samples appearing at the output of first sum-and-dump 510 at the sampling rate 
ftfMl are passed to the 2nd order recursive filter section. This section is made up of 

latches or multi-bit flip-flops 524 and 532: gain units 5 12, 518, 520, 528 534 and 
536; adders 516, 522, 526, 530, /and 538; and subtract 514. The input samples 
arriving from the first sum-and-dump are coupled to gain units 512 and 536 The 
15 recursive section operates at the rate f r /Ml, as determined by the clock signal at 
frequency f r /Ml applied to latches 524 and 532. The filtered output is taken from the 
output of adder 538, which is coupled to the input of the second sum-and-dump 
section at adder 542. The second sum-and-dump section, formed by adder 54- 
accumulator latch 544, and output latch 546, operates in a manner identical to that of 
20 the first sum-and-dump section, except that the input sampling rate is f r /Ml, and the 
decimation parameter is now M2. The output of latch 546 represents the final 
Iowpass filtered and decimated result, which exhibits a sampling rate of f r Ml/M2. 

The structure of the 2nd order recursive section has been chosen so that *ain 
units 512, 518, 520, 528, 534, and 536 may be implemented without us~in» 
25. multipliers. The structure is of a type especially suited for narrowband Iowpass 
filters wherein the gain coefficients may be coarsely quantized without causins 
excessive frequency response distortion. Such a structure is discussed, for 
example, in the article, New Recursive Digital Filter Structures Having Very Low 
Sensitivity and Roundoff Noise, IEEE TRANSACTIONS ON CIRCUITS AND 
3D SYSTEMS, Vol. CAS-22, No. 12, pp. 921-927 Dec. .1975. In the filter of the 
present invention the process of coefficient quantization has been taken to the limit, 
i-e., the coefficients have been quantized to "one bit". Thus each gain coefficient is, 
in fact, a power of two. Those of ordinary skill in the art will realize that to 
implement a power-of-2 gain in binary arithmetic requires only a simple bit-shifting 
35 operation. In practice, this bit-shifting can be implemented by an appropriate routing 
of data lines. Thus, no hardware per se is heeded to implement the gain units. This 
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facilitates the overall filter implementation, as it is reduced to simply a combination 
of adders and latches. 

The values of the gain coefficients and the decimation parameters Ml and 
M2 should be chosen to yield good SNR performance and an output sampling rate 
as close as possible to the ideal low value of 15 kHz. As an example, in one 
embodiment the design parameters are as shown below in Table 1. 



ro: 
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Table 1 



Where ML = 1 12, M2 = 5 



With an input sampling rate fj equal to 9.548 MHz, as provided by the 
20 quadrature mixer/decimator circuit of Figure 4b with f s = 38.192 MHz, the final 
output sampling rate is 17.05 kHz, which is quite close to the minimum 15 kHz rate. 
The z-transform of the transfer function of the second order recursive section is 
given by: 



25: 



H(z) = 



7 2 - 157.^1 



z 2 -l.5625z-.875 



30 



35 



where l/z corresponds to a delay of 1/85,25 kHz. 

The magnitude vs. frequency response of the overall filter is shown in 
Figure 6c. Compared to the response of the simple sum-and-dump filter it exhibits a 
more rectangular aspect with improved stopband rejection. Accordingly, the SNR 
performance is improved considerably. By virtue of first decimation by Ml = 1 12. 
the sampling rate of the recursive and second sum-and-dump sections is a desirably 
low 85.25 kHz. 
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* ° f 16 f is of * e mb o d i m ;„r ~ 

m^unator 322 described above, assuming that A/D convener 106 is a 1" 



bfc**. ^ "J^tum edtothearchitectureandoperationsofDSPliz a 
<*k diagram of which is shown in Figure 7. In the preferred embodiment DSP 
112 mcludes, for each of four parallel channels, a search processor and 
10 ^ng/data processor. The search and tracking/data p„c«Z^ 
^coupled tomecorrespondingchannelprocessmgunitofDSPUO. ForexaTp.e 

dsp no p ™:«rrr t ~ coupied to ^ 1 p ~s -*» * 

fifLdl , , * 40X161 1 Pr ° mpt PI «*■ PQ ' ^m lowpass 

IS ^ 7 ^ 324 ' are eaCh C ° UpIed t0 SearCh Passer 702 and to 

15 belong/data processor 704. The channel 1 differential signals DI and DQ, output 
from lowpass fdter/decimators. 326 and 328, are each coupled to tracking/data 
processor 704. The channel 1 C/A code control signal, output from trackin /d a a 
processor 704, « coupled to C/A code generator 330. Similar connections are made 
between channel 2 processors 712 and 714 and processing unit 340 of DSP no- 
20 between channel 3 processors 722 and 724 and processing unit 360 of DSP i 10- 
and between channel 4 processors 732 and 734 and processing unit 380 of DSP 
110. All eight search and tracking/data processors are coupled to control 
microprocessor 1 14 via the bi-directional data bus 750. 

Figure 8 shows in detail the basic operations performed by tracking/data 
25 processors 704, 714, 724. and 734. Briefly, each processor provides Costas loop 
earner recovery and carrier demodulation (Doppler removal), delav-locked loop 
processing for C/A code delay control, automatic gain control (AGC), and data bit 
timing recovery and data detection. An input/output (LO) bus interface is also 
included to facilitate communication with microprocessor 1 14. 
30 Costas loop carrier recovery and carrier demodulation are performed by the 

loop comprising complex mixer 802, lowpass filters 808 and 810, mixer 842, 
carrier loop lowpass filter 844, and quadrature variable frequency oscillator (VFO) 
846. Prompt I and Q signals, PI and PQ .arriving from DSP 1 10 are coupled to one 
input port of complex mixer 802. The other input port of complex mixer 802 is 
35 coupled to the output of a quadrature variable frequency oscillator (VFO) 846. which 
supplies sampled cosine and sine waveforms at the negative of the recovered carrier 
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(Le., Dopplcr) frequency. These signals are denoted cos 0 and sin 0, and appear at 
outputs 846a and 846b, respectively, of quadrature VFO 846. In effect, complex 
mixer 802 performs a complex multiplication of the two "complex" signals PI + jPQ 
and cos 0 '+ jsin 0. The product resulting from this multiplication comprises two 

05 output signals: an in-phase ("real") component {PI cos 0 - PQ sin 0}, and a 
quadrature ("imaginary") component {PI sin 0 + PQ cos 0}. These signals are 
applied to AGC amplifiers 804 and 806, respectively, which serve to scale the 
signals to a relatively constant average power level, thereby mitigating the effects of 
signal fading. The derivation of the AGG control signal will be described shortly. 

10 The outputs of AGC amps 804 and 806 are coupled to lowpass filters (LPF r s) 808 
and 8 10, which act to improve the loop signal-to-noise ratio (SNR). The outputs of 
LPFs 808 and 8 10 are multiplied together by mixer 842, and the product is applied 
to carrier loop LPF 844. The output of loop filter 844 is the recovered carrier 
frequency, fooppler* which is coupled to the frequency control input 846c of 

15 quadrature VFO 846. 

The control loop thus formed may be viewed as a Costas loop, or as a 
complex- valued squaring loop operating at zero frequency. Such loops are well 
understood by those skilled in the art. One Costas loop which would function 
satisfactorily in conjunction with the GPS receiver of the present invention is shown 

2o and described in a text by Holmes, Coherent Spread Spectrum Systems, Wiley 
Interscience 1982, pp. 121-207. In the locked condition, the "complex" sinusoidal 
output of quadrature VFO 846 closely follows the incoming prompt signal in phase 
and frequency (actually the negative of said phase and frequency). Thus the output 
of complex mixer 802 represents the fully demodulated data signal, centered 

25 substantially at zero Hertz, with zero phase angle. In practice, the bandwidths of the 
identical "arm" filters 808 and 810 are adjusted for the best loop SNR performance 
consistent with the dynamic tracking requirements of the loop. Carrier loop filter 
844 is also designed with tracking requirements in mind. It preferably includes an 
integrator so that the phase may be tracked with zero error for any Doppler 

30 frequency. Quadrature VFO 846 may be conveniendy implemented using a standard 
ROM lookup table technique, described, for example, in the article "A Digital 
Frequency Synthesizer", by J. Tierney, C. Rader, and B. Gold, in IEEE 
Transactions on Audio and Electroacoustics, Vol. AU-19, No. 1, March 1971, pp. 
48-56. 

35 When the carrier recovery loop is locked, the in-phase demodulated signal 

appearing at the output of amplifier 804 is the baseband 50 bit/sec data signal ready 
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for detect™. This signal is coupled to bit timing recovery block 892 and data 
detector 890. The output of bit timing recovery unit 892 is a bit timing sign* 
synchromzed tothc transition instants of the received data signal. This biftimkg 
signal 1S applied to data detector 890, instructing the detector as to the optimum 
tuning of b« demons. The output of data detector 890 is the recovered sateZ 

8T 8 ^r aSC *-T Ch " 5601 10 miCrOPrOCCSSOr W 1 14 via I/O bus interf I 
882. ™*fc<a™ngmfomation^ 

m navigation calculations. The workings of bit timing recovery unit 892 and Z 
de^ctor 890 will be well understood by those skilled in the ^ZZl lt^ 
^ d ******* "f*» « orally described in a a text by Holmes, CoZ2 
Spread Spectrum Systems, Wiley Interscience 1982, pp. 564-620 

The differential signals arriving from DSP 1 10 are processed in a manner 
surular to that performed on the prompt signals, in order to provide C/A code delay 
locking. In-phase and quadrature differential signals DI and DQ are input to one 
portof complex mixer 822. The other input port of the mixer is provided the cosine 
and sine outputs of quadrature VFO 846. Complex mixer 822 operates in a manner 
identical to prompt complex mixer 802, and produces in-phase and quadrature 
output signals (DIcos 0 - DQsin 0} and {DIsin 0 + DQcos *}. The in-phase and 
quadrature outputs of complex mixer 822 are scaled by AGC amps 824 and 8^6 
and then lowpass filtered by filters 828 and 830. Filters 828 and 830 are identical to 
prompt path filters 808 and 8 10. The in-phase prompt and differential signals output 
from filters 808 and 828, respectively,. are multiplied together by mixer 872. The 
corresponding quadrature signals output from filters 8 10 and 830 are multiplied by 
mixer 874. The products of mixers 872 and 874 are summed by adder 876' whose 
output is coupled to delay-locked loop (DLL) lowpass filter 878. The signal output 
from DLL filter 878 represents the receiver-derived estimate of relative C/A code 
delay, which is sent to nP 1 14 for navigation calculations, via bus interface 882 and 
bus 750. The code delay signal is also sent via switch 880 (position 0) to the code 
control input of the appropriate C/A code generator of DSP 1 10. 

The delay-locked loop thus formed by the tracking processor of Figure 8 
operating in conjunction with DSP 1 10 is known as a modified noncoherent code 
tracking loop. Such a loop is described, for example, in the article "A Modified PN' 
Code Tracking Loop: Its Performance Analysis and Comparative Evaluation by R. 
Yost and R. Boyd, in IEEE Transactions of Communications, Vol. COM-30, No. 
5, May 1982, pp. 1027-36. For the present invention, the loop described in the 
article has been adapted to operate at zero frequency in a complex-number (I-Q) 
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fashion. The loop dynamic response is controlled by DLL filter 878, which 
preferably includes an integrator so that the delay may be tracked with zero 
steady-state error for any Doppler shift 

In order to maintain, at a constant level, the dynamic characteristics of the 
05 various receiver tracking loops (carrier, delay, and bit timing), the tracking 
processor of present invention is provided with AGC. An AGC control signal is 
derived using squarers 852 and 854, adder 856, subtracter 858, and AGC lowpass 
filter 860. The in-phase and quadrature output signals of filters 808 and 810 are 
applied to squarers 852 and 854, respectively. The squares of the two signals are 
10 calculated, and then summed by adder 856. The output of adder 856 represents the 
instantaneous power of the I-Q signal output from LPFs 808 and 810. A constant 
Kp, corresponding to the desired average power level, is added to the negative of the 
sum output by adder 856, by subtracter 858. The output of subtracter 858 
represents a power "eiTor" signal, which is filtered by AGC lowpass filter 860. 
15 Filter 860 preferably contains a perfect integrator. The output of AGC filter 860 is 
the AGC control signal, which is coupled to AGC amps 804, 806, 824, and 826. 
.The AGC loop thus formed acts to maintain, at the constant level Kp, the average 
power of the post- arm filter I-Q signal. The bandwidth of AGC filter 860 
determines the dynamic response of the AGC loop, and should be chosen in 
20 accordance with the characteristics of expected signal strength variations. 

The above described operations pertain to the tracking/data processor as 
operated in a "track" mode. During the search procedure (to be described shortly) 
which precedes tracking, the tracking/data processor is configured to operate in a 
"search" mode- This is accomplished by jjP 1 14, which sets switch 880 to position 
25 I- With switch 880 in position 1, the tracking/data processor simply acts to relay 
code control information input from pP 1 14 (via bus 750 and bus interface 882) to 
the C/A code generator of DSP 1 10. This information includes the code number 
(i.e., which of the 18 satellite codes is to be generated) and the C/A code delay, 
which is incremented, or swept, during search. In the search mode, tracking and 
30 data detection operations may essentially be idled Upon completion of search, ^iP 
114 sets switch 880 to position 0, and sends an initial Doppler frequency estimate 
(obtained during search) to quadrature VFO 846 via input port 846d. The tracking 
process then begins. 

The tracking/data detection operations shown in Figure 8 may be 
35 conveniently implemented using off-the-shelf, programmable DSP ICs. The 
possibility of such an implementation is enhanced when the sampling rates of the 
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signals input from DSP 110 are minimised. A sampling rate of 17 OS kH, f u 
Ft PQ, D!. and DQ s.gnais U sufflciendy tow to Iw each of m f tt ^"^ 

05 one^ Tp IC iS the,iPD7720 * NEC (NEC Electronic^ U s A 

05- One Natick Executive Park, Natick, Mass. 01760): JA ' 

Modifications or additions to the tracking/data detection operation, shown in 
Ftgure 8 may occur to those skilled in the art Such additions ntigh. Z^Lt 
detector, for the various loops, or the outputting of additional info^o" 

!0 • l-mZ* ' T* 5 "" S "* ^ fre " A "*« *X <>fa 
rT^r" ,emenaa0niSrtaI * e ^ ss ^°P«^may b«=asU and 
fiextbly modified to accommodate the requirements of different receivl appti^s 

detect' to m thc back Sround discussion, before signal tracking and data 

derecnon may begtn, each satellite signal must be searched for and acquired 

13 sTTp ."I 9 ° m " ma Pr0b,em ° f fOT • OA cod* 

^-Ftgure 9 a shows the cross correlation between a receive OA signal and a 
locally generated OA code, as a function o, OA code delay. This is equivalent to 
dte re.anve magmtude of a despread (prompt) data signal as code de.ay is varied 
TOe cross correlation essentially exhibits a triangular shape two chip periods (about 

20 i T T Z * base- which appeare once every C/A code < 1023 chi ^ °< 

20 1 mtlhsec). The desired operating point is at the peak of the triang.e, where 
despread signal power is maximum. NormaUy, upon receiver startup , he correct 
code delay is uncertain, i.e., the location of the correlation triangle is unknown. To 
find the correct code delay, the receiver typically sweeps or increments through all 
poss.bte code delays until a high despread signal level i s detected. De.av trials ate 
25 usually spaced about a half chip period (.5 usee) apart, so that the correction peak is 
not substantially missed. 

Ftgure 9b illustrates a second dimension of the search problem, the frequency 
uncertainty. When correctly despread, the C/A coded signal is collapsed down to 
simply a 50 bit/sec data signal. The carrier frequency of this data signal is uncertain 
30. due to the Doppler shift produced by satellite motion and receiver clock drift In the 
present receiver, which operates at essentially zero frequency, the carrier frequency 
of the despread (prompt) signal may vary between ^7.5 kHz and -7.5' kHz. 
Filtering, provided by the lowpass filter/dec imators of DSP 110, is used to remove 
noise outside this band. However, the signal-to-noise ratio within this bandwidth is 
35 still not sufficient to allow for reliable detection of the signal, nor to allow for 
locking of the tracking loops. Accordingly, to improve the SNR, further filtering 
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must be performed which breaks the Doppler frequency range into several narrower 
bands, each of which is checked for signal presence. 

Figure 9c summarizes graphically the two-dimensional nature of the search 
process. The desired signal may be viewed as lying at some point within a 

05 two-dimensional search space, whose dimensions are time (code delay) and 
frequency (Doppler). To locate the signal, the space is divided into many smaller 
sub-units, or bins, of dimensions At by Af, each of which is checked for presence 
of signaL Typically, to allow for reliable detection, At is about .5 usee, and Af is on 
the order of i kHz or less. Since the total dimensions of the space are 1 msec by 15 

10: kHz, the total number of search bins is on the order of 30,000. 

Most prior receivers perform search in a sequential manner, examining each 
bin in the- search space one-at-a-time. The usual procedure is for the receiver to 
make a guess as to the correct code delay, setting the C/A code generator 
accordingly. Then, for that code delay, all frequency bins are examined for presence 

15: of signaL If no signal is found, the code generator is incremented to the next code 
delay and the process is repeated. The detection method typically employed in 
checking for signal presence can be loosely described as a, variable- frequency 
average power "meter" followed by a threshold decision. For a given code delay/ 
Doppler frequency trial, the C/A signal is multiplied by the locally generated C/A 

20 code, and mixed with a sinusoid whose frequency is equal to the current Doppler 
guess. The resultant signal is filtered by a filter whose bandwidth is commensurate 
with the frequency search interval Af. If the desired signal is present in the bin being 
searched, the output of said filter will contain a despread 50 bit/sec data signal. The 
output of the filter is squared and averaged to obtain the average power, which is 

25 then compared to a threshold. If the average power is greater than the threshold, the 
signal is deemed present in the current bin. 

Otherwise, the signal is deemed absent, and the search proceeds to the next bin. A 
search technique which typifies a conventional practice is described in a paper by 
Spilker, J.J., Jr., Global Positioning System; GPS Signal Structure and 
30 Performance Characteristics, THE INSTITUTE OF NAVIGATION, Vol. 1, pp. 29-54, 
1980. 

Due to the large number of bins to be searched (around 30,000), 
conventional receivers require many minutes to find a signal if a sequential search 
, method is employed. Additional time may be incurred if the identities (i.e., codes) 
35 of the satellites within view of the receiving antenna are unknown. Methods to 
shorten the search time have been devised but are quite expensive to implement. 
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GeneraUy, these are based on the idea of examining several code delays 
sunultaneously for each Doppler frequency in order to speed up the process For 
example, one technique employs multiple despreading circuits, each with its own 
code generator and average power detector, operating in parallel. 
05 ^ The receiver of the present invention employs an improved search 
technique, which has been developed in conjunction with the overall receiver 
^hnecture. Briefly, the technique involves the use of a Fast Fourier Transfer 
(FFT) aided average power detector, which allows all Doppler frequency bins to be 
_ examined simultaneously for each code delay trial, thereby drastically reducing 
10 search ume. Most or all of the required processing may be conveniently 
implemented using programmable DSP IC's. 

In one embodiment, the FFT-aided search procedure is implemented using 
search processors 702, 712, 722, and 732 of DSP 1 12 operating in conjunction with 
microprocessor 1 14. In this embodiment, overall search control is provided by iiP 
15 114, and the bulk of the processing, including FFT calculations, is performed by the 
search processors. 

In summary, the improved search method works as follows. For a e ive n 
code delay trial, search processor 702 (or equivalently, processor 712 72"» or°73?) 
inputs, upon command from ^P 114, blocks of N prompt sample pairs (PI PQ) 
20 from DSP 110. For each block of N I-Q sample pairs, a complex-discrete Fourier 
transform (DFT) is performed using an X-point FFT. In effect the FFT 
accomplishes complex mixing of the I-Q sample stream by M complex sinusoids 
whose frequencies are equally spaced across the Doppler range, followed bv 
lowpass filtering and decimation of each product signal. The N mixing and filtering 
25 operations are performed simultaneously, in a highly efficient manner, bv the FFT. 
In practice, the filtering is equivalent to an N-sample sum-and-dump. Decimation bv 
N is thus performed inherently by the FFT on each of the filtered product sisnals. 
Therefore, for each block of N I-Q input samples, there are N complex, or" I-Q, 
output samples, each corresponding to a different frequency bin. The power of 
30 these N FFT outputs is then found and averaged. The M powers are found by 
forming the squared-magnitude of each complex output (i.e., the sum of the squares 
of the real and imaginary, or I and Q, components). The M power quantities are 
then averaged over M FFTs. 

At the conclusion of the average" power measurement, processor 702 
3 5 determines the largest of the N power quantities, and passes this value, along with 
the corresponding frequency bin number, to ^P 114. Microprocessor 1 14 compares 
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this maximum power to a threshold Depending on the result of this comparison, 
search either ends or proceeds to the next code delay, as in conventional search 
procedures. 

In effect, the improved search method described above implements N 

05 average power detectors, operating in parallel, for N frequencies spread evenly 
across the entire range of possible Doppler frequencies. Since N frequency bins are 
checked simultaneously for each code delay trial, the required search time is reduced 
by the factor N compared to conventional sequential search procedures. The 
frequency resolution Af obtained is equal to the sampling rate of the I-Q sample pairs 

Xa in P ut fr° m DSP 1 divided by NV If the input sampling rate is denoted f^, then 
the frequency' resolution is Af=fjq/N, and the frequency range covered is from 
-fiq/2 to -rfiq/2. The FFT parameter N should be chosen to provide adequate 
frequency resolution, and to allow for convenient FFT implementation (since FFT 
complexity, e.g., memory requirements, depends on N). The averaging parameter 
M and the threshold value affect various detector characteristics, such as the 
pro babili ties of detection and noise falsing, as is - well understood in the art. 

A flow diagram, outlining a suitable procedure for the. maximum power 
measurement performed by search processor 702, is shown in Figure 10. The 
procedure may be advantageously implemented using a programmable DSP IC such 

7q as the |iPD7720 manufactured by NEC, as mentioned above. In this embodiment, 
the FFT parameter N is equal to 32. For an input sampling rate of 17.05 kHz, the 
frequency resolution Af is 533 Hz, and the range of frequencies covered is from 
-8.525KHz to +8.525KHz. 

The maximum power measurement is initiated at item 1002 upon command 

25 from jiP 1 14 (via bus 750). The first step, carried out in item, or sub-routine, 1004. 
is to clear the 32 power averaging accumulators ZSUM(k), k=0 to 31, to zero. 
Next, an input sub-routine 1006, 32 consecutive prompt sample pairs (PI, PQ) are 
input from DSP 110 and stored in memory. These I-Q sample pairs may be 
regarded as complex numbers of the form: 



30 



x(n)-PI(n)+jPQ(n) 
where n is the sample time number ranging from I to 32. 



35 
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The 32-point complex FFT of the input sample sequence is com «,., a ■ 
sub-outine 1008. Mathematically, the complex values X(k k 7 T " 
calculated, where C ».*W. *-0 to 31, are 

05 32 

X(k) = 2 X (n) exp{j2rcnk/32} 
n=l 

10 ^VT 16 k CO T POndS t0 freqUenCy * aCC ° rding » the rel ~P fluency - 
10 k-Af. For example, k=0 corresponds to zero frequency.. Negative fre que nclel 
correspond to values of k between 16 and 31, since the spectrum of a sanded 7Z 

lorirT CXamPle ' COrTeSP ° ndS » thC ^ ^vel ^ 

algonthrns are available for computing the X(k). An "in-place" algorithm* 

15 '"IT T ^ ^ ^ ° f » is -U known in the ar. One 

X5 such algotithm is described in a text by Oppenheim and Schafer, DigUa^ 
Processing, Prentice Hall 1975, pp 284-336. 

_ Once the FFT is compu ted, the squared-magnitude, or power, of each of the 

T'Z T 15 f ° Und ^ SUb — ™ ^er quantities, denied 

Z(k), are found by squaring and summing the real and imaginary parts of the Xfk) 
20 M^ernatically, Z(k, , (Real X(k) } 2 . {Imag X(k)} 2. £ ™ 

1012, the power quantities Z(k) are summed, or accumulated, into the power 
accumulator registers ZSUM(k). Each ZSUM(k) is updated by adding Z(k) to the 
present value ofZSUM(k) to obtain the new value. 

The routines 1006, 1008. 10 10. and 1012 are executed M times, as provided 
25 by deczston block 1014. Thus. M sets of 32 prompt sample pairs are input and 
transformed, and M sets of powers are computed and summed The numerical 
quantities present in the ZSUM accumulators at the conclusion of this process 
represent the 32 averaged power outputs corresponding - to the 32 frequencv bins 
The largest of the 32 ZSUM values, denoted ZMAX. is determined bv routine 1016 
30 along w lt h the corresponding value of k, denoted KMAX. KM AX is that value of k 
for which ZSUM(k) = ZMAX, i.e., ZSUM(KMAX) = ZMAX The values ZMAX 
and KMAX are returned to nP 114 (via bus 750), and the maximum power 
. measurement procedure of search processor 702 is terminated. The processor then 
waits for command from nP 1 14 to begin another measurement. It may be noted . 
3 5 that, since the sampling rate f iq is necessarily greater than the Doppler frequency 
range, several values of the frequency variable k may correspond to frequencies 
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outside the range of possible Doppler frequencies. In the example illustrated here, 
with f iq equal to 17.05 kHz the FFT frequency range is -8,525 kHz to +8.525 kHz, 
whereas the Doppler range is -7.5 kHz to +7.5 kHz. Therefore in practice, the 
power calculations pertaining to those frequencies outside the Doppler range may be 
05 skipped, and the maximum average power selection routine 1016 may ignore those 
"extraneous" frequency bins as well, since there is no possibility of signal presence 
there. The frequency search range may be similarly restricted if there is some a 
priori knowledge of Doppler frequency. 

A flow diagram of the search control procedure executed by microprocessor 
Iff 114 is shown in Figure 11. The search routine is entered at block 1102. An 
initialization process is carried out according to routine 1104. First, the appropriate 
tracking processor (704 in this channel 1 example) is set to search mode by sending 
a signal causing switch 880 to be set to position 1. Next, the C/A code number 
corresponding to the desired satellite is sent through track processor 704 (and bus 
E£ 750) to the appropriate code generator of DSP 110. Finally, an initial C/A code 
delay value may also be sent to the code generator, if desired. 

Upon completion of the initialization procedure 1104, jiP 114 issues a 
command to search processor 702 to perform a maximum average power 
measurement (item 1106)- Microprocessor 114 then waits (item 1108) for 
ZO measurement results from processor 702, namely, the maximum power value 
ZMAX and the corresponding frequency bin number KMAX. In decision block 
1 1 10, ZMAX is compared to a threshold. If ZMAX is less than the threshold, the 
search at the present code delay is deemed unsuccessful. In such case, the code 
delay is incremented (block 1 1 12) by sending an incremented code delay to the C A 
25* code generator via track processor 704, and the search measurement is repeated. 

The loop comprising blocks 1106, 1108, 1110, and .1112 is executed 
repeatedly until the quantity ZMAX is greater than or equal to the threshold. When 
_ this occurs, search is deemed successful, and the estimated Doppler frequency 

KMAX«Af is transferred (item 1 1 14) to track processor 704 (specifically to port 
30 846d of quadrature VFO 846). Track processor 704 is then set by item 1116 to the 
track mode, by setting switch 880 to position 0, and tracking commences. The 
search procedure executed by microprocessor 114 is then terminated (block 1 1 18). 
until such time as it is reactivated by some higher level routine. Modifications to the 
improved search procedure described above may occur to those skilled in the art. 
35 The particular partitioning of functions between search processor 702 and \xP 1 14 
might be changed, for example, by having the search processor perform more of the 
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search control functions. As with the tracking/data processor, an important 
advantage of the programmable implementation made possible by the receiver 
a^tecture disclosed herein is that search operations may be easily and fle xib y 
modified to accommodate varying receiver requirements. 

05 , , * TZ*; ^ f0reS ° ing **** ^ apparatUS «* me *<* ^ receiving 
a pluralxty of GPS OA coded signals, wherein a substantial part of the proce 

may be advantageously implemented using digital circuitry, J*SZ 

proved method for rapid initial acquisition of signals is employed. While the 

10 ^ atUrcS ° f inVenti ° n beCn SWn * ^ <* Nation, other 
10 modulations and changes will occur to those skilled in the art. It is, therefore to be 
understood that the present claims are intended to cover all such modifications and 
changes as fall within me true spirit and scope of the present invention 
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1; In a spread- spectrum receiver, for receiving one or more 
direct-sequence-coded signals, each of which may be Doppicr shifted from its 
nominal carrier frequency, wherein for each signal a code delay/Doppler frequency 
search procedure must be carried out prior to activation of delay-tracking and 
05 carrier-recovery loops, an improved method of Doppler frequency search for each 
desired signal, and for each code delay search trial,, comprising the steps of: 

a. inputting N consecutive digitized samples of a signal equal to the 
filtered product of the down-converted received signal multiplied by the 

Id locally-generated code of the desired signal, the code delay of said 
locally-generated code being that of the current code delay search trial, 
said filtering passing, at least, the band of frequencies corresponding to 
the entire range of possible Doppier-shifted carrier frequencies of said 
down-converted received signal; 

15 

b. calculating the N-point discrete Fourier transform of said block 
of N samples to produce N complex-valued transformed output samples, 
each corresponding to one of N equally-spaced search frequencies in the 
band from zero Hertz to fs (or, equivalently, in the band from -f s /2 to 

2o +f s /2), f s being the sampling frequency of the samples input in step a; 

c. calculating the square of the magnitude of each of said N 
complex- valued output samples to produce N power quantities, 
corresponding to the N said search frequencies: and 

25 

d. performing steps a through c M times, summing the M power 
quantities produced for each of the N said search frequencies, to produce 
N quantities proportional to average power, the maximum of said average 
power quantities exhibiting indicia of the presence and approximate 

30 doppler frequency shift of said desired signal. 
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2. 



In a spread-spectrum receiver, for receiving one or more 
du ect-sequence^oded signals, each of which may be Doppler shifted fromT* 
nonunal earner frequency, wherein for each signal a code delay/Doppler fr'enct 
-arch procedure must be carried out prior to .activation of dela^tracSg "d 
camer-recovery loops, an ^improved method of Doppler frequency search Z u 
desiredsignal^dfor each code delay s^^^^^^ 

r, JLTf 1 * N COnSetUtiVe ***** ^Ples of a signal equal to the 
filtered product of the down-converted received signal multiplied by the 
locally-generated code of the desired signal, the code delay of said 
locally-generated code being that of the current code delay search trial 
sard filtering passing, at least, the band of frequencies corresponding to 
the entire range of possible Doppler-shifted carrier frequencies of said .. 
down-converted received signal; 

b. calculating the N-point discrete Fourier transform of said block 
of N samples to produce N complex-valued transformed output samples 
each corresponding to one of N equally-spaced search frequencies in the 
band from zero Hertz to fs (or, equivalent^, in the band from -f s /2 to 
+f s /2), f s being the sampling frequency of the samples input in step a; 

c calculating the square of the magnitude of each of said N 
complex-valued output samples to produce X power quantities, 
corresponding to the N said search frequencies: 

d. performing steps a through c M times, summing the M power 
quantities produced for each of the N said search frequencies, to produce 
N quantities proportional to average power, 

30 e. selecting the maximum of the N average power quantities, and 

calculating the corresponding frequency; and 
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f. comparing the maximum average power quantity to a threshold, 
the fact of said threshold being exceeded signifying that the current code 
delay is approximately correct (local code approximately time-aligned to 
received code), with the frequency of maximum average power found in 
05 step e indicating the approximate Doppler frequency shift, and the fact of 
said threshold not being exceeded signifying that the current code delay is 
not approximately correct, and that another code delay should be tried. 
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X In a GPS receiver, for receiving one or more C/A coded signals each 
of wtach may be Doppler shifted from its nominal carrier frequency T^Tf 
^^-odedelay/Dopplerfrequencysearch 
pnor to activation of delay-trackmgafldcarrier-recovLi^ c 

mal, comprising the steps of: cm ' search 

r.u l™!^" 8 N C ° nsec,,tive »»Pl« of a signal equal to the 

ocaUy-generated code of the desired signal, the code delay of said 

locals-generated code being that of the current code delay search rial 
satd filtering pas sin S , at least, the bartd of frequencies c orresponding to 
tfte enore range of possible Doppler-shifted carrier frequencies of said 
down-converted received signal; 

15. 

b. calculating the N-point discrete Fourier transform of said block 
of N samples to produce N complex- valued transformed output samples 
^ each corresponding to one of N equally-spaced search frequencies in the 
band from zero Hertz to fs (or, equivalent^, in the band from -f s /2 to 
f s being the sampling frequency of the samples input in step a: 

c calculating the square of the magnitude of each of said N 
complex-valued output samples to produce N power quantities 
corresponding to the N said search frequencies; and 

25. 

d. performing steps a through c M times, summing the M power 
quantities produced for each of the N said search frequencies, to produce 
N quantities proportional to average power, the maximum of said average 
power quantities exhibiting indicia of the presence and approximate 
30 Appier frequency shift of said desired signal. 
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4. In a GPS receiver, for receiving one or more OA coded signals, each of which 
may be Doppler shifted from its nominal carrier frequency, wherein for each signal a 
code delay/Doppler frequency search procedure must be carried out prior to 
activation of delay-tracking and carrier-recovery loops, an improved method of 
05 Doppler frequency search for each desired signal, and for each code delay search 
trial, comprising the steps of: 

a. inputting N consecutive digitized samples of a signal equal to the 
filtered product of the down-converted received signal multiplied by the 

Id locally-generated code of the desired signal, the code delay of said 
locally-generated code being that of the current code delay search trial, 
said filtering passing, at least, the band of frequencies corresponding to 
the entire range of possible Doppler-shifted carrier frequencies of said 
down-converted received signal; 

15 

b. calculating the N-point discrete Fourier transform of said block 
of N samples to produce N complex- valued transformed output samples, 
each corresponding to one of N equally-spaced search frequencies in the 
band from zero Hertz to fs (or, equivalently, in the band from -f,-/2 to 

20 f s bein g *e sampling frequency of the samples input in step a; 

c calculating the square of the magnitude of each of said N r 
complex-valued output samples to produce N power quantities, 
corresponding to the N said search frequencies; 

25 

d. performing steps a through c M times, summing the M power 
quantities produced for each of the N said search frequencies, to produce 
N quantities proportional to average power; 

30 e. selecting the maximum of the N average power quantities, and 

calculating the corresponding frequency ; and 
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.5. An improved substantiaUy digital receive ^apparatus for the reception of one or 
more dtrect-sequence-coded spread spectrum signals,- which exhibit substantially 
equal earner frequencies, except for Doppler shifting, said apparatus comprising: 

0 5 (a) first means for coupling a wideband radio frequency (RF) signal 

contaming said spread spectrum signals to said receiving apparatus; 

(b) second means, including preselector and down^onverting means 
coupled to said first means, for selecting a composite signal containing said 

10 spread spectrum signals from said wide-band RF signal and for 

frequency-translating said composite signal to an intermediate frequency; 

(c) digitizing means, coupled to said second means, for convening 
said translated composite signal to a digital composite signal, at a 

15 predetermined sampling rate, wherein the predetermined sampling rate is 

determined according to the relationship: f s = 4 f { wherein f; is said 
intermediate frequency and f s is the predetemuned sampling rate; 

(d) first digital signal processing means, coupled to said digitizing 
20 means, for separating and despreading each of said spread spectrunvsignals 

from said digitial composite signal, and for deriving a prompt and auxiliary 
late^arly signal corresponding to each spread-spectrum signal: and 

(ft) second digital signal processing means, coupled to said first digital 
25 signal processing means, providing for each desired signal, code delay 

search and tracking and carrier search and tracking processing for each of 
said doppler shifted signals, and data recovery processing. 
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6. The apparatus as recited in claim 5 wherein said second digital signal 
processing means includes means for Fast Fourier Transform (FFT) aided 
spectral power measurement 

' .... ■ ■ 

05 7 ♦ The apparatus as recited in claim 5 further including microprocessor means 
for receiver control functions. 
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8 ; An improved method for digitally receiving one or more 
direct-sequence-coded spread spectrum signals, which exhibit substantially 

^ uenc - es ' except for simn8 > said method co « 

' (a) coupling a wideband radio frequency (RF) signal containing said 

direct-sequence-coded spread spectrum signals to said receiving apparatus; 

s <b) selecting a composite signal containing said spread spectrum 
signals from said wide-band RF signal; 

-■ ■ ■ ■ 

(c) converting said translated composite signal to a digital composite 
signal, at a predetermined sampling rate; 

(d) separating and despreading each of said direct-sequence-coded 
spread spectrum signals from said digital composite signal, including 
generating a prompt pseudorandom code signal and a late/early 
pseudorandom code signal and deriving a prompt and auxiliary late-early 
signal corresponding to each direct-sequence-coded spread-spectrum signal- 
and ' 

(e) code delay searching and tracking and carrier searching and 
tracking the doppler-shifted signal correpsonding to each desired spread 
spectrum signal, and recovering a data signal. 
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9- The method as recited in claim 8 wherein said separating and despreading st< 
includes the steps of: 



055 



(a) generating a prompt pseudorandom code signal and a late/early 
( pseudorandom code corresponding to each desired spread-spectrum signal; 

(b) multiplying the prompt pseudorandom code by the digital 
composite signal to generate a first digital product signal; 

10: (c) filtering the first digital product signal to generate a prompt 

signal; 

(d) multiplying the late/early pseudorandom code by the digital 
composite signal to generate a second digital product signal; 

is: 

(e) filtering the second digital product signal to generate an auxiliary 
late/early signal; and 

(f) outputting said prompt and said auxiliary late/early signals for 
20- further processing. 



25 



30 



35 



WO 87/01540 



PCT/US86/01765 



41 



05 



10 



10 . The method as recited in claim 8 wherein said code delav searrh , 
said doppler shifted carrier frequency search and track ~ ^ ' 

desired signal, the steps of: P mcIudes ' for « 

(a) estimating a code-delay for the desired direct-sequence-coded 
spread spectrum signal, and generating a pseudorandoj ] Z ^ 
corresponding to said delay; wgnai 

_ ft) measuring the spectral power of the prompt signal for the current 
code delay generated by said separating and despreading step: 



(c) comparing the power values measured in said measuring step to a 

SETT ■"T 01 * — Wherein 3 ^ *~ ««^'ng said 
hreshold value mdtca.es the proper code deiay and doppler shifted carrier 
frequency; and 

(d) incrementing the code delay and repeating steps a xhrough d until a 
proper code delay and doppler shifted carrier frequency is found! 
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11 . An improved substantially digital receiver apparatus for the reception of one or mo 
GPS C/A coded spread spectrum signals, which exhibit substantiaUy equal carri 
frequencies (except for Doppler shifting); said apparatus comprising: 

0 5 (a) first means for coupling a wideband radio frequency (RF) signal 

containing said GPS C/A coded spread spectrum signals to said receiving 
apparatus; 

(b) second means, including preselector and down-converting means. 
La coupled to said first means, for selecting a composite signal containing said 

GPS C/A coded spread spectrum signals from said wide-band RF signal 
and for frequency-translating said composite signal to an intermediate 
frequency; 

15 (c) digitizing means, coupled to said second means, for convening 

said translated composite signal to a digital composite signal, at a 
predetermined sampling rate, wherein the predetermined sampling race is 
determined according to the relationship: f 5 = 4 f ^ wherein fj is said 
intermediate frequency and f s is the predetermined sampling rate; 

20 

(d) first digital signal processing means, coupled to said digitizing 
means, for separating and despreading each of said GPS C/A coded spread 
spectrum signals from said digitiai composite signal, and for deriving 
prompt and auxiliary late-early signal corresponding to each 

25 spread-spectrum signal; and 

(e) second digital signal processing means, coupled to said first digital 
signal processing means, providing for each desired signal, carrier search 
and tracking and code delay search and tracking processing in response to 

3 0 said prompt and auxiliary late-early signals, and data recovery processing. 
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12. The apparatus as recited in claim 11 
wherein said second digital signal processing 
means includes means for Fast Fourier Transform 
(FFT) aided spectral power measurement. 

13. The apparatus as recited in claim 11 
further including microprocessor means for 
GOS navigation calculations. 
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1 4 - A method for the reception of one or more GPS O A coded spread 
spectrum signals, which exhibit substantially equal carrier frequencies 
(exceptforI>jppiershifting);saidmethod comprising the steps of: 

05 (a) coupling a wideband radio frequency (RF) signal containing said 

GPS OA coded spread spectrum signals to said receiving apparatus; 

(b) selecting a composite signal containing said GPS C/A coded 
spread spectrum signals from said wide-band RF signal and for 

10 frequency-translating said composite signal to an intermediate frequency; 

(c) converting said translated composite signal to a digital composite 
signal, at a predetermined sampling rate, wherein the predetermined 
sampling rate is determined according to the relationship: f s = 4 f { * . erein f- 

15 is said intermediate frequency and f s is the predetermined sampling rate: 

(d) separating and despreading each of said GPS spread spectrum 
signals from said digitial composite signal, including generating a prompt 

C/A code signal and a late/early C/A code signal and deriving a prompt and 
20 auxiliary late-early signal corresponding to each GPS spread-spectrum 

signal; and 

(e) code delay searching and tracking and carrier searching and 
tracking the doppler-shifter signal correpsonding to each- desired GPS 

25 spread spectrum signal, and recovering a data signal. 
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15. The method as recited in claim 13 «,h~.,« ^ 
^spreading step includes thestepsof "* and 



(a) generating a prompt C/A code 

• signal; 

0» multipiying the prompt C/A code by the dfetei 
enerateafirstdi^^..-..,. . ^ 0y me ^S 1 ^ comppsi 



o-— ««« a prompt t/A code Siffnal a h j — » , 
OS "nespond^^ 



to generate a fust digital product signal; composite signal 

10 ^ fiIl ering the first dieitai 

signal; ^ Pr0duct S1 * nai » generate, a prompt 



(d) multiplying the late/early C/A cexV h« *u j 
signal to generate a second digital 

product signal- composite 

:r onddi8iaipro<ittct " gn "" o ~ --^ 
20 ^;r=:;s said;promp,Mdsaid -- i - 
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16. The method as recited in claim 13 wherein said 
code delay search and track, and said doppler shifted 
carrier frequency search and track process step includes 

for each desired signal, the steps of: ' 

* 

05 (a) estimating a code-delay for the 

desired GPS spread spectrum signal, and 
generating a C/A code signal corresponding 
to said delay; 

10 (b) measuring the spectral power of 

the prompt signal for the current code 
delay generated by said separating and 
despreading step; 

15 (c) comparing the power values measured 

in said measuring step to a predetermined 
threshold value, wherein a power value 
exceeding said threshold value indicates 
the proper code delay and doppler shifted 

20 carrier frequency; and 

(d) incrementing the code delay and 
repeating steps a through d until a proper 
code delay and doppler shifted carrier 
25 frequency is found. 
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